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EFFECT OF COMPEE^IBILITy AT HIGH SUBSONIC VELOCITIES 
ON THE MOMENT ACTING ON AN ELLIPTIC CYLINDER 
By Carl Kaplan 


SUMMARY 


An extended form of the Ackeret iteration process is utilized 
to calculate the compress ihle flow at hi^ subsonic velocities 
past an elliptic cylinder. The angle of attack with respect to 
the direction of the undisturbed stream is assumed small and the' 
circulation is fixed by the condition that tho trailing end of the 
major axis.be a stagnation point. The expression for the mcament 
acting on the elliptic cylinder is derived and shows a first-step’ 
improvement of the Prandtl-Glauert approximation. In addition^ a 
second-step iinprovement is obtained in the Prandtl-Glauert approxi- 
mation for the lifting force, doting on 'thd' elliptic cylinder. By 
means of .these two results it' is possible to calculate ■fche effect 
of compressibility on the position of the center of pressure as ■ 
a function of tho thickness coefficient and of 'the stream Mach 
number. Tables and corresponding graphs are included to illustrate 
numerically the theoretical resiilts derived. For example, it is 
found that, for an elliptic profile of thickness coefficient 0.15 
and stream Mach number 0.8o, the center of pressure moves rearward 
a distance 2.6 percent of the chord from its position in the 
Inccmpreseible flow. 

. . INTRODUCTION 


The present paper is concerned mainly with the calculation of 
the effect of compressibility at hi^ subsonic velocities on the 
moment acting on an elliptic cylinder. The method used is an 
Iteration procedvu*e, credited to Ackeret, which proceeds from the 
Prandtl-Glauert approximation as the first step and successively 
improves i-fc in a systematic manner. . The details of the Ackeret 
iteration process have been described in reference’ 1 and, therefore,- 
only material essential to the present paper will be repeated. 

The main purpose of the Ackeret iteration me'thod Is to 
linearize the nonlinear partial differential equation (for the 
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Telocity potOT,tial or the stream function) that Governs the steady 
tvfo-dlmensional flow of a p^fect ccmpresei'bl© fluid. This 
linearization is aacongjllsh^ hy assuming the development of -Uie 
stream function say, to he of" the foma 


f=-m + + ^l/2(z^y) + - (i) 


idiere U is the velocity of the undisturhod stream and X and Y 
are the rectangular Cartesian coordinates of the physical flow 
plane. Eqiiation (1) is essentially a development of the stream 
function around a uniform stream in the negative direction of^tho 
X-axis. For the purpose of defining or controlling the iteration 
procedure, the function is regarded as smlT ccairpared with 

the preceding function ■ 4 ^. and the derivatives have a similar ^ 
relationship. Then the total index decides the order of tiio term; 
for ’epca^le, .'ij/g is ,of the same order as or 4 / 24 /g* l‘he 

accuracy of this iteration method clearly depends .on -the degree to 
which the assumptions are satisfied. In the case of slender "bodlea 
without stagnation points, the first few steps may he ospected-to 
yield a good result . In the case of hodios with stagnation points^ 
the accuracy of the calcxxlaticcta obviously depends on the number 
of terms derived, each new term reducing the extent of the 

region of inaccuracy in the nsi^borhood of the stagnation point. 

Ip the; treatment of the. vhrious. eq.uationS that result, frcm . 
the llheardzation; of the fuhdnmentai differential . equation’ by' 
means of the Acieret Iteration process, it is ppnyeniont tp iptyrp- . 
duce an affine transfoimatioh of the ' plijuicaJ. flow plane.. This 
affine transfcamiatlan reduces the differential oquaulona bo bo 
solved to a Laplace equation and to Poisson equations. In the 
performance of -this slnpllflcatipn,. the statement of the boundary 
condition at tho solid %y moans ’ of the volocity potential bocesmos 
very complicatod. Fortunately, however, tho statomont of-^bho 
boundary condition by moans of -toe stream function, namely, 4r = 0 
at the- solid, • is invariant for tho af f ino , transfofrmation; thoro- 
fexTO, the use of tho stfoam function throughoilt tho ‘analysis of 
the present paper is to be prefbfrod. Tlio choico ‘of . the ollipso 
as tho bolid' boundary is dictated "by. tho ’pTopo3;*ty .i^t an af fine... 
distortion of ah ollipso leads to another oU^sa; thorofore., the . 
analysis can bo conduoted entirely. in tho affinely distorted plane 
•and the results thus obtained 'lihiod -to the actiial olliptlc profilo 
by means of slmplo corrospondonce relations. 
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MOMEOT K 3 EMJLA 


Specifically, tjie problem treated herein is to obtain an 
improvement of the Prandtl-Crlauert approximation of the effect of 
ccanpresslbility on the moment acting on an elliptic cylinder set 
at a small angle of attack in a uniform stream. Let Z denote 
the physical flow plane, z the affinely transformed plane, 
and z ‘ the plane of the circle Into ■which the affinely distorted 
profile is mapped by a conformal transformation. (See fig. 1 .) 

As in the calculation of the restiltant lifting force given in 
reference 1 , it is a great labor sa'ving device to choose a large 
circle in the z ' -plane to ccxrrespond to the con-trol contour in the 
physical Z-plane during the calcxilatlon of the moment and also to 

choose as independent variables the polar coordinates 
of the z '-plane, with 


-it 

. , . ' . . 2 ' = Ee ^ . (2) 

where ^ = £+ iri and R is the radius of the conformal circle. 
(See fig. 1(0).) 

Since the large circle in ihe z ' -plane corresponds to a large 
control ellipse in the physical flow plane Z, the eapression for 
the moment must contain, in addition to the usual momentum integral, 
a term involving the Integration of the pressures around 'the control 
ellipse. This additional term is necessary beca'ttse the normal 
vector to an ellipse does not pass throu^ its center. Q?ho general 
vector expression for the moment in a ccmpressiblo fluid with 

respect to the origin, obtained frcm reference 2 , is 



where brackets and parentheses denote vector and scalar products, 
respectively, and 

r radius vector from, origin 

n unit normal vector ' 

4 velocity vector of fluid 
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h 


ds element of length along control contour 

p pressure of fluid 

p density of fluid 

The positive direction of the unit normal vector n is from the 
control contour toward the origin, and the line integrals are taken 
positively counterclockwise around 1316 control contour in the 
physical flow plane. The unit tangent vector t and the unit 
normal vector n thus form a ri^t-hand frame; hence, a positive 
value for the mcment corresponds to a counterclockwise rotation 
(fig. 1(a)). 

It is easy to verify that equation (3) can he rewritten in 
the form 


M. 



- lv)% dZ + 



, 2 
or 




where u and v are the congjonents of the velocity vector along 
the X-axls and T‘'axls^ respectively, a^ 


Z = Z + iY 


= :ZZ: + Y^- 


2 2 2 
g = u + V 


Note that if the fluid is inc'ampresslhle 


1 2 ' 
p + a Constant 


according to Bernoulli ’s eguation; therefore, the second integral 
vanishes identically and yields the usual Blasius formula for the 
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moment . 
becomes 


By use. of the adiabatic relation — 



equation (^4-) 


M, 





:(u - 



+ 


1 fl/PlJ 

2 




(5) 


■vjhere the subscript 1 refers to the starting conditions at infinity 


Telocity of undisturbed fliild at infinity 
Mach number of •undisturbed stream at infinity 

velocity of sound in undisturbed fluid 

ratio of specific heats at constant pressure and constant 
Toliune, for air 7 = 1.4 

For -fche puirpose of calcvilating the line integrals indica-fced 
in eq.uation ( 5 )^ it is necessary to express the integrands as 
functions of the independent variables of the z ’-plane • 

In "the case of the elliptic profile in the z -plane, the confomal 
transformation to a circle of radius E ■with center at the origin 
is 


and 

U 


7 


z = c cos(t + IX) 


(6) 


■where c is the semifocal distance and X is defined by any one 
of -tiie follo-wlng geometric characteristics of the ellipse: 
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a ss o coeh \ 


"b = o sinh X 


E = ~ce^ 

2 

■where a, h, and B are, respectively, semlmajor axis of ellipse, 
eemiminor axis of ellipse, and radius of confonnal circle. Now, 
the affine •bransfonnation used in connection wi'fch the Ackeret 
iteration process is 


X = X 


Y « fxy 


•where 




It follo'ws -tliat 



and 


7 _ 
Z = 



dZ = 


1 




■where a "bar indicates conjugate-octn^lex quantities. Since ■the 
control contopr is a large circle in ■the z '-plane, n = Constant 
and d^ = d^ = dg. Then, hy use of equations (2) and (6) and the 


relation 


E = 


the expressions for Z and dZ . on the control 


ellipse hecome 
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z = 




Tl+\ 




n) f ir + “ 


2t]+2X, 


(7a) 


and 


dZ = - 


cs 


TJ+X 


(1 + n): 


2 'e‘ 


2ti+2X 




2 ■ Q^r\+'^ J \ z 


f2X j 


dz 


(7t) 


■where z ' = 
Similarly, 


2 1 - 11 ^ 2 1 - H ^-2 1 + - 

r* = 7,7. ^ j 2- + 'g + ^ZZ 

h h 2 


r = ZZ 


and on the control ellipse, t] = Constant, 


2 

dr 


- 1) -^'^y ■ ^) (i^2 ■ ' 


4t)+4X 


(8) 


Expressions for 


P , 


iv) 


and Pj^/P as functions of the 


variable z '(= e~^^), expanded in powers of l/e’^, are given by 
equations (46) and (47), respectively, of reference 1. By the use 
of these eq.uations, together with e<iuations (7) and (8), it is easy 
to evaluate the ri^t-hand meitiber of equation (5) by noting that 
only terms involving dz'/z' contribute to the line integrals. 

The result thus obtained is 
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„„ = - pV - 1 ) - 1 )(- } 

•where a Is ■fch.e angle of attack in the affinely distorted plane z 

and o = (7 + l)(n^ - l}* If the quantities a, "b, c, and- a 
are replaced by a'^ h, c*, and a* of the actual profile in 
the Z -plane according to the correspondence equations of refer- 
ence 1, naaiely. 


a a a ' 


1 , 
h a —b * 


2 .2 . ti - 1 ,2 

e SB c ■+ • ■• r— b 

\r 


(9) 


1 . 

a = -a' 


the moment about the origin on the actual ellipse becames 



Me = rtp^U^a'c'^lx 


i«p a’b '^(a + 4)— 

8 1 V- 


- 1 


( 10 ) 


Now, for an IncompreBsible fluid. 



or ■ . 

p = 1 • 

Mj^ = ct*c 


. and 



NACA TN No . 12l8 


9 


Therefore, 





( 11 ) 


Tdiere t' Is the thickness coefficient h'/a' of “Uie actual 
elliptic profile in the physical flow plane. 

Eq,uation (11) represents a first-step lagjrovement of the 
Prandtl-Glauert approximation and reduces to that result in- the 
limiting case t* — >-0. This improvement, however, is inccanplete, 
for, as can he observed from eq.uaticn ( 10 ), the second term on the 
ri^t-hand side is of the third order (that is, proportional 

to a'V^) tar» of «iat order ere caatrituted ty ^, 3 . 

Since only the first two terms and ‘W®!’® 9-erived in refer- 
ence 1 , it is necessary to determine the third term in order 

to obtain the complete first-step ijnprovement of the Prandtl-G-lauert 
approximation for 

EETERMETATION OP 

In order to obtain the third term expansion for the 

stream function 


xj/ = - Tjy + x}/^ ■*"^2 ■*■••• (12) 


it is first necessasry to obted.n the expression for inclusive 

of third order terms, in the nei^borhood of the undistxnrbed stream. 
Analogoxis to equation (19) of reference 1, this expression is 



X . - fr - - *)(^ - *)■[<> , >. . 0 . «f - 


s 




^ •*'^Y 
xf- -2[J^ 


-|(ii^ - i)^ |i5 + (7 + iV^[(27 + 15) + 3(7 + - 1)]} 


gy 
J tj 3 


(13) 


Vh^ th^ ejpresaltftja for ij/ arid p. fe givan "by eqxiatiQna".;(l2) and (I3) ere snbatl'tiited into the 
■basic differential equation [■ . 


5!^V+ —f— 0 

^xVp ■ hY\p In 


ilk) 


nrxPi teiBa of the third order In the dorlTiatlToe of are collected, the following differential 
Bqiiatian. for ie obtained; 
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r 


’•^3X1 + 


V = 4^ - . [i(, . !)(,= - 1)= - h1^' 


+15 

D® 


?1XX 


(^ •,'^r, r i\/2 Al’''32'*^2XX’*'’*'l5X+2y ^/o AP / ^2 A1 

- - ly ^2 + (7 + 1) - l)\ 2(|i 2 - ljji+ (y + 1)^,1'^ - ijj 


1 +ix+iY+ur 

TJ^ 



This differential eq,uatlon. can be expressed In a convenient fcam for solution by making use of the 
affine transfannation 


X s> I 


y « iy 

and by introducing a nev stream function ^jr*, ■where 


ijf = (jUijf* 
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Then equation (15) ’becoHea 


H 

ro 


'^ 3 n * 'l^ 3 JT= * (•‘^ ’ ^ *|<^ * «l*V ' 


J..^ _ t . r~ . n \ /^. -^ _ T ll AffS- %\/K‘ * ^IrJt ^lic. \ r\/ S -i\ri r^_ ^ \ / Q. ■ .. 


"‘l) 1^57 +7) - (7 + I)(n^ - 1)] - if |i5 + (r + Dli^ [(27 + 15) + 3(7 + 1) (n^ - 1)] 


V 


”[7 + 4) + (7 + l)(ix^ " l)j b ■*■ ■*■ ‘ ^)J1 


(16) 


Again, aa in reference 1 , it •Hill be found that the mathematical analyals -will be considerably 

einplifled by workinK ■with a ncmana lytic cotanl.ex potential ■w*^(z,z') instead of its imaginary 

o‘ ' ■ " ' 

part As shown in inference 1 , is the Imaginaiy part of an analytic functicn -w^^^Cz), 


VhnT'nnfl 


is 'the insginary part of a nDnanalytic function 


XJ^ ( IT ■z^ . 
- gx-7--- 


T+. T»mfrh "hn pinmhnn'f v.na 


that "the real parts of ■these ccai^ilei functions are notl to be interpreted as velocity potentials but 
only as functiaos that render ■the analysis elegant and sinple- The foUorwing identities can be 
easily verified; 
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■ ^15) 

^7 = ^ ^iz) 

^l2y “ 2 (”izz ^ ^Izz) 

= fe(^lzz " hzi) 

- % - "^2z) 

= I«P*(^2z ^ ^2z) 

%j = \(^2z ■ '"25 - ^2z + ^2z) 

= I-P-i(^22 " ^25) 

^23y 2V^2 zz ” '*2 zz " '"2zz ’”2zz) 

= I -P -1(^222 - -Waiz) 

'i'2s3: = il(""2zz + ^^z5 + ^2zz ' "^2zz " ^^2zz " ^25z) 
= I.P.(v2zz + 2 v2zz + W2 zz) 
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where the aster iBk has been dropped. ^Phen 


='■^("12 " ""iz) (^Izz " "^Izz) 

• %r) 

'I'lx'fejy + 'I'lsy'l'ai: * ■ "'SSe) 

■* ^ + "2S) 

■*■ 'h.xx^2y ” ■*■ ^lz)^*^'('^2zz ^'^Szz ^£zz) 

+ i^-^-^lzz ■ \zz)("£z ■ "Sz) 
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By means of these Identities and "by substituting for frcja 

equation ( 29 ) of reference 1, equation ( 16 ) can he "written as 


“ ^22) (^izz - hzz) + 2(^2z + ^2z)(^lzz + %ii) 

|X i 

-(,! + £) + \z) *H'^zz - '*'2m)(" 12 - %a) 

+ ^ L(jl^+ 1) + (cj + 3)(^^ - 

+ “ + 1) + jc + (c + 2) (a + "Wj2 ^izz 

+ 1 [(cr +- 2 )%- - 1) + ct(ct + 3 ) (li^ + l)]-^i2^i2-«-i22 ^ 

where 1® imaginary part of w^. I^om equation (33) of 

reference 1 



(ff + 4 )w3_w^ 


+ 



■where F(z) is a .function, "the form of "tdiich is decided hy the 
boundary conditions but need not be given explicitly at this point. 
When the expressions for "the derivatives of "Wg are inserted into 

equation (l7)j then 







H 

o\ 




ct((7 


8 


- 1) 


Izzz 


C‘^'l)(*l%)j’'aja+5 |s(<r + l)+ff%® + l)+ jff+(i + 4)(|r + J^)J(^®-l)|•VJJ^w. 


i£i!^/' 2 
8 


12 


i|. 4 kj („2 - l) + „ (|»s + l)j 


Izz 


(18) 


This dlfforsiitial ©quation can "be integrated by tDepoctloa 'Wlthonit difficulty. Thus the general 
solution, is 
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* I 






cr + 4/ 2 




8 


+ ff(tf + 3)^ 2.-. 


\^U " - l) “ i) 

+ cr (^^ + l)]zwi 2 ^ - ^ |8((t + 2) ^ + + 2( 3 o + 8) (tr + £)] - 1 )| 


2 , 

vij, dz 


+ a^(z) + Gg(z) 


(19) 


viisre Ctq^(z) and Gg(z) are orbitorary analytic functions, raapectively, of only t and z 

to be dfitemilned by the boundanTr conditions. The boundary conditlona to be satisfied are that 
at the surface of the ellipse, t] « 0, 




(20a) 
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and at infinity, = », 


^ S\lr 

— = .r-i =. 0 20(b) 

5x 3y 


and that the trailing end of the major axis be a stagnation point. 
It may be noted that the result represented by e< 3 .mtion (19) is 
not restricted to an elliptic profile but is valid for an arbitrary 
solid boundary. 

Again, as in reference 1, the most direct way to impose the 

boundary condition, s I.P.w = 0 at the surface of the solid, 

3 3 

is to utilize the "polar" variable £ of the z' -plane. Thus, for 
the elliptic profile according to equation (39) of reference 1, 


= c cos (t + 11) - 2E cos £ - 2Ha (sin £ + £) 


> i-^1) 


vr 


Iz 


= V 


l£ dz 


= 1 - 


2R 


sin t 


2R 


cos 


t + 1 


C sin (^ + il) C Bln (£ + il) 


•where, from equation (6) 


z 


C COB (£+11) 
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Ju6t as in reference 1, it is a sli^ile matter to supply the 
functions of L needed to satisfy the boundary condition at the 
surface of the elii^e in the 2 -plane. For example, at the surface 
vhere ri = 0. cos (T " il) = cos (^ - iX) . By use of eq.uations (21) 
eq.uation (I9) then becomes 

[=os^(t - 1 ^) |- oos2(£ - 

H* 


+ - l)c [cos (f- iX) - cos - iX.)]vjL2^’2, 


+ — ^cos (f - iX) jcoB (f^ - iX) - ^ cos t 

- (sin f + f)j - cos (^ - iX) |cos - iX) - 2^ cos ^ 

- 2|a (sin t + 

+ ^ |^(c + 3) + (3cr + 8) (\i^ - l)j jcos (f - iX) - cos I 


- (sin f + f ) - cos (^ - iX) cos ^ ^ 

a + hf 2 ,\2R sin t 2R cos t + 1 
+ ^ - 1 ; a . = — - — 

8 [_c Bin (t ” iX) ° sin (^ - iX) 

- M t ^ ^ cos ^ + 1 I j, 

° sin (^ - iX) o sin (^ - iX)_ 


(eq.uation continued on next page) 
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+ “1)0 cos(f^-lX)-— cos f-'^a (sin f+t) 

o L ® ^ 


" cos - iX) +— cos ^ +^a (sin t + t) 
c c 


F- 


^ - " ({1^ - l)c^ |pos (C - iX) “ cos f - (sin t ^'E)] 

cos (t-lX)-^coB t "—a (sin ^ 


, r ^ - _ r ^ ^ ^ _.g igL .... 

16 I ■ L ° sin (t - iX) 


2E cos t + 1 
- -— ^ — -. — „ ■ „ , 

° sin (f- IX) 

2R^ cos t 1 
c sin ( ^ - 1>.) 


1 


cos (^ - iX) 


1 + 


^ t 

c sin U - IX) 




” ^ gin. ^ 

L ° sin (I - iX) 


+ ^cc _£2L_Li_i_ + ^ sin I 2 R^ 


cos 


t + 1 


° sin ( ^ - iX) ® sin (t ~ iX) c sin - iX) 




crc 

48 


+ ( 7 <J + 8) |cos (f - iX) - cos ( ^ - iX)^Wj^2^ 


k. fp 

16 


^8(cf + 2)2 + [ofS + 2(a + 2) (30 + 8)J ((a“ - 


2E sin 


L 


° Bin (t - IX) 


2 R cos t + 1 2E sin f 
- — a 2 

c sin (^- iX) c sin (^ - iX) 


2R cos t + 1 

— — — -- — I I 

c sin ( t - iX) J 


t 


US.P. 


2 <iz 


H 


( 22 ) 



uhere 


^ I 




(is.p. 


cos (^ + iX) + 1 


+ ^ cos (t - i\) +'~a-|2 Bin (^ •’ IX) +2^ coeh X -- slub 2X log 


4E‘ 


2 ' 


. c*- 


c^ 


t 


cos (t + IX) - 1 
COB {£ + IX) + 1 


“ ei siiili X log Bln (J + IX) 


}..{c 


Iffi/ 2Si . ,2, , coah X + 1 

cosli X + — (rtCt" 1) alnh X leg 

c 2 cosh X 1 


+ ooeh X “ ^^-tt|2i elnh X + 2jt coah X - ~ Blnh 2X log ^ ^ ^ 

^2 [_ 2 COBhX-1 


21 slrth X log ("1 Binh X) 


]) 


and S.P.. the lowor limit of integration, denotes the stagnatlm point. By meana of the follows 
fotnnula (eetifltion ( 44 ) of reference 1 ) for the complex yelocity (vlth rejprd to ^3 only) In the 


physical Z-plane; 


-73) = ■ *3) 


(24) 


ro 

H 
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equation (22) can "be showi to fulfill the "boundary conditions, not 
only at the solid surface, tj = 0, "but also at Infinity, ti = 
moreover, the trailing end of the major ajdLs (g= jc, q = 0) is a 
stagnation point. In order to satisfy the boundary conditions, 
however, a number of singularities of the nature of doublets have 
been introduced into the field of flow. Those umrelcome singu" 
larltles are caused by the factor l/sin ( ^ ~ il) and are located 

at the exterior points z' = iEe^ or = il and ^2 “ ^ 

They are removed by the addition of doublets in such a manner that 
the sum of the residues at a pole is zero. The images, moreover, 
of these superimposed singularities in the conformal circle of 
radius R must be Included in order to insure that the boundary 
conditions are preserved. As an example, consider the expression 


H(0 

sin •* iX,) 


( 25 ) 


where the function H(t) is regular everywhere in the fiiilte region 
exterior to the circle of radius R (t^ = 0)* Then, in .order to cancel 
the residues at the poles a iX, and - n +iX. and to preserve 

the boundary conditions at the solid and at infinity, the following 
expression mvist be added to the rl^t-hand side of equation (22): 


- ^(iX) ^ot cot 


- II 


By means of this expression the additloml terms can bo easily 
obtained in order that w^ be regular everywhere in the finite 

region exterior to the conformal circle of radius R in the z'*plano 
or to the elliptic profiles in the z-plane and Z-plano. An 
examination of equation (22) immediately yields the following equa- 
tion for H(t); 
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H(D = 


- l)e^ [^- sin ^ + a (cos t + 1)] [i(^ ) + “ cos 
- (a+k) J + ^ ^(ff + 2 ) 2 + 2 (ff + 2 ) ( 30 + 8 )](n^ - l)j^ [- sin 

in 


+ a (cos ^ + l)Je^ 


ilS.F. 


2 *2 j. 


"IZ 


( 2 ?) 


^ere, fi-om equation (^3) of reference 1, 


I'(t) = ooB (C - iX)v.^ - 2(0 + k)ic Bihh X Bin - (o- + 4)w3^w^j 
2 


+ 2lD(cf + 4)a (sin t+S) 


( 28 ) 


Now, the introduction of the foregoing singnlarities induces 
a finite velocity at the trailing end of the major axis of the 
elliptic houndary. This Velocity, obtained hy moans of e<iuation (24) 
and ezprossion (26) , is 




1 

4ic sinh X 


<[s{iX,) +H(-iX)3 


1 

.2 X 
cosh — 


+ 



It + iX) + H(?i ” H> 


)] 



(29) 


Then, by means of equations (2?) and (28), with terms of hi^er 
order than the third neglected 
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5if. - H.„ ■ > • ''V - >) 


2/ 2 ^ b 

a + "b 


la 

128 


(3 - log 4 


)(n^ - 1) [8(0 + zy 
[ff® + 2(0 + 2) (30 + 8)J(n^ - l)| 


a + b 


In order to maintain jt, r\ = O) a stagnation point, the 

folloving expression, the imaginary part of ■which vanishes 
for T) = 0 and for r) = oa, 

0 (cr + 4 ) 2 (h 2 - 1 ) - ^(3 - 

• 1 . [of^ + 2(a + 2)(3cr 
must he added to the ri^t“hand side of equation (22) • Finally 

Q 

then, the complete expression for w, is given hy •fche right- 

- 1 

hand side of eq'uatifm (22) and the expressions (26) and ( 30 ). 



EQUATIONS FOE 



AND 



The ocmponents ^l 2 and v^ of -fche velocity of tlie com- 

pressihle fl'uid in the physical flow plane Z are calciaated hy 
means of equation (24). Since, for. the purpose of the present 
paper, calculations are performed along an ultimately infinitely 
large contour, the developments for u^ and v^ in "bhe nel^hor- 

hood of infinity are sufficient. Thus, h’y ■use of the ccairplete 
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expression for the development of the ccairplez velocity 

— - — ''T, I in the nei^horhood of infinity is given hy 

+ ia6 + Bg log -“|-~ 

PlUi^J \i 3J \a + Dy z>eTi 2 ay\^a + hy ^ .2^^ 


- ioc/^-Y ^+iaD|-^Y + E Kct/™Y 

\f + t/ q2ti \a + h/ 2 >^e^n f. \a + h/ 


\?+V iz.Vt 


(31) 


■where 

z' = e-i^ ' • • • • ■ 

and 

A == |-(o + - 1 )^ +—(3 - log h){y,^ - 1 ) |8(cf + 2)2 

+ [a^ + 2(0 + 2) (3a + 8)J(Y " 

- 1 ) [^(ca^ + 7a + If-) + - 1 ) (l3cf^ + + 32 )! 

Bg = "^4^^ ■*’^) + 2 )'=^ + jff^ + 2(a + 2) (3a + 8)J(^^ - i)j> 

**'',** 

C = + ^0(o + 8)(Y " 1) 
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+ k)(/ - if 

E = ^(cr + 4)^(|a^ - l)^ 

Erom eq.uation ( 31 )^ it follows easily that 
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2 ? 


Then 





i^om eq.uation (Ito) of reference 1 and equation (33) # therefore, the 
development in the nei^torhood of infinity of the complete con^lex 
velocity, inclusive of berms of the third order, is given hy 
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Prom eqiiation (3^)> 



Now, by definition 


Ei. - 

Pj_U UT 


and 


p^U U 2 
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Hence 



n-nri "by means of eq.uatlon (17) of reference 1 for ~ and eq[ua 
tion (35), 



CALCULAIION OF THE LIBT 


In a ccoqprossi'ble flow as In an inconipressilJle flow, the lift 
is given by 
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where P is the circulation round the profile and where, hy 
c 

definition. 


r = ax + V dY) = e.p. 1 Q)(u - iv) az 
c 




Wow, dZ is given hy equation (7h), and from equations (34) and (36) 


i(u - iv) . -1 - lKa+ ^ |l + - 1) (c + 1.)^ 


Then, because only terms that -involve dz’/®' contribute to the 
line integral, it follows that 


r = 

c 


2 

4id?Uti a 


1 + 


1 / 2 



or if A is replaced by its definition (see equation (31))^ it 
follows that 


= 4jtRUti^a 1 + - 1) (c + 4)| + H- 4)^(ta^ - 1)^ 

+ |(3 - log 4) (n^ - 1 ) ^(CT + 2)2 + [ct 2 +-2(a + 2) (3 <t + 8)] - l) 


e2_ 

(37) 


If, according to the correspondence equations (9)# t, a, and E 

are replaced by and E' + ', respectively, then 

U Z]i 

for the actxial ellipse in the physical flow plane. 
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= In® TJa V + 2:SJa’b ' |(l - n) - l)(cr + If) + ^) “ 0 

+ |(3 - los ii-) (n^ - l) {8(«J + 2)^ + [f^ + 2 ((t + 2) (3ff + 8)] (ti^ - 
Since tbe circulation in the case of incourpreBsibl© flow is 


= kaR'Ua* 


the ratio TqIV^ or- Lo/% given by 



■where t ’ ie the thickness coefficient b ’/a * of •fche actual 
elliptic profile in •fahe physical flow plane. Equation (38) 
represents a secorui-step lmpro-7ement of the Prand-bl-G-lauert 
approximation and reduces to "that result idien t* —*>0. In refer" 
ence 1 a first-step ing)rovemsnt of the Prandtl-G-lauert approxi- 
mation for -fche ratio 'was calculated and is represented by 

•fche first t'wo tearms on the ri^t-hand side of equation (38) • 

Table I sho-ws values of the ratio for the first-step and 

second-step improvements, for various values of the thickness 
coefficient t' and the stream Mach nvmiber (■with 7 = l.k 

for air). Eigure 2 shovTS the corresponding graphs with as 

abscissa and Iq/Lj as ordinate. An examination of -fchese graphs 
sho^ws that below ■the critical s^bream Mach number M^y ■fche main 
effect of compressibility is already given by ■fche Prandtl-G-lauert 
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term and the first-step improvanent and that larj^e differences 
between the first -step and second-step irtpravements do not appear 
xtntil well above the critical stream Mach number. 


In reference 2 the ratio ’L^I'L^ was calculated for an elliptic 

cylinder by the method of Poggi. This result, restated in the nota- 
tion of the present paper, is 



1 + 


K 


2il + 2t' 


log 


2 - t 


,2 


1 - t’ 1 - t' “ 1 + t* (1 . t')' 



i 


1 + t 
1 - t 


T 


\ 



log 


\ll -h t' + - t’ 

fT7‘ - 


+ 


(39) 


and must agree with eq.uation (38) insofar as tlie terms common to 
the two developments are concerned.- If^ than, eqiiation (38) is 
expanded according to powers of Mq^ and eq.uation (39) is expanded 

according to powers of t', the two expansions are found-to agree 
and yield 


= 1 + ^2^ + + |(1 - log 4)M^^t '^ + . . 


CALCULATION OP TEE MOMEIKT 


Por the purpose of calculating the moment, the following two 
eq.uations are needed. Prom eouations (7) 
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and from eq.iiations (3^) ( 36 ) 


P 






la 


+ (1“r)' 


^]- 


la ..I 15 


z'“j e^n 

.i(,^-l)(a..^)(^y][l- 


2 . .2 


p) z 


+ (1 + ji) 


■M 


+ . . . 


By means of these eq.uations and eq,uations (8), (35); and { 36 ) 
equation (5) for the moment about the origin yields the 

foUoving result: 


M. 




+ 1)Bq_ + klB 2 log 



(1 + H)Z'^ 
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AgaiHj replacing a, "b, c, and a according to the correspond- 
ence equations (9) yields 


M, 


. = jtp£U^ a'c - ^gpjU^a' ^ — - '(a + 4)b a* --—- 



where and B2 have been replaced by their definitions. (See 

equation 31*) Now, for an incon5>resBible fluid, 




The ratio for the actual elliptic profile in the i^iysical 

flow plane therefore becomes 



32 + (1^^ - 1) +;2)^] - 



+ - 1) |5^ + 2(cr + 2) (3or + 8)j y log 


1- t 


,2 


(40) 


Equation (4o) represents- the complete firstrstep improvement 
of the Prandtl-Glauert approximation for the ratio of mcanents 

and reduces to that result in the limiting case t ' . — ^ 0 . Again, 
as in the case of the lift, the ratio was calculated for 

an elliptic cylinder- by the method of Poggi (reference 2). This 
result, restated in the- notation ■ of the present paper, is 


M, 


= 1 


+ •^1'" 
2 ^ 


,2 


(1 - t‘) 




T 1 

log-, 



(41) 
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JuBt as in the case of the lift, eq.uations (40) and (41) must agree 
insofar as the terms common to the two results are concerned. Thus, 
if eq.ua tion (4o) is eapanded according to powers of and equa- 

tion (4l) is estpanded according to powers of t', the escpansions 
are found to agree and yield 

Me . h, 2 ovt 2+ ,2 ^ M 2. ,2 1 ^ 

= X 4* ** ^ ^X ^ Xog 4* • • • 

Table H shows values of the ratio M^/M^ calculated by means 

of equation (40) for various values of the tl^ickness coefficient t' 
and the stream ifech nmber and figure 2 shows the corresponding 

graphs with as abscissa and as ordlmte. 

Contrary to the lift, Tdiich is a localized vector, the moment, 
is a nonlooallzed vector, the loagnitude of ■which depends on the 
point about 'vhlch it is taken. In the present paper, this point 
is the origin of coordinates. If, now, the moment about -the 
origin 0 is denoted. by Mqq, the moment about any other point P 

in the plane of flow is then given by 



where r is the length of the perpendicular dropped frcaa the 
origin 0 to the line of action of the lift vector "Ihrou^ 

the point P. If this espresslon for the moment is examined in 
relation to the moment about the same point P in an incompressible 
fluid, it will be seen -that the ratio of moments again 

begins 'With the Prandtl-G-lauert approximation but "tha-t -the hl^er 
terms of -bho approximation depend on the point P about 'vhich 
the moments are taken. Pigure 2, consequently, should not be 
used to compare the Various moment curves ■with the Prahdtl-G-lauert 
approximatlonj rather, the significant result is •the compressibility 
effect on the movement of the center of pressure from its position 
in an inccaj^ressible fluid - a quantity -that is independent of the 
point about which moments are taken. • 
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EB^11CT..0F COMERE^IBILirY ON POSITION OF CENTER OP PRESSUEE 


If - C^' and 0^’ denote the distance of the center of pressure 

from the center of the ellipse in the canipressihle and IhcoD^ressihle 
fluids, respectively, tlien 


£c ^c ^1 
C, == Mi 


Also, Ci =? “(1 " 


therefore. 


Qc ’ Qj ^ 1 - t H _ \ 

2 a k ^ J 


By means of this formula and equations ( 38 ) and (40) it is possible 
to calculate the effect of ccmrpresBihility on the position of the 
center of pressure for various thiclmess coefficients and stream 

Cc - Ci 

Mach number’s. Table HI shovs values of the ratio , the 

2 a 

negative values indicating movement toward the center of the elliptic 
profile • Figure 2 . shows the corresponding graphs- with the stream 


Mach number . M^^ as abscissa and the ratio — -r 

1 2 a 


in percent 


chord as ordinate. , Note that in each case at. some high subsonic 
stream Mach number, the movement of the center of pressure reverses'. 
(See table III where sign changes from negative to positive.) This 
peculiar behavior of the center of pressure is probably caused by 


the term log in the eq.uation for the moment M^, and indicates 
the need for additional terms in the expansion for the stream 
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function xjf to insure greater accuracy in the range of hi^ sub- 
sonic stream Mach numbers. 


Langley Memorial Aeronautical Laboratory 

national Advisory Committee for Aeronautics 
Langley IFleld, Ta., October 2k, 1946 
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«1 


(^‘Irst-step linproTemeiit) 

— 1 

L^IL^ (secand-Btep iEprovemficit) j 

O.JjO 

.20 

•30 

-1*0 

.45 

•50 

.60 

.65 

•70 

•75 

.80 

•90 

1.0050 

1.0206 

1.0483 

1.0911 

1.1198 

1.1547 

i.1974 

1.2500 

1.3159 

1.4003 

1 - 5 U 9 

1.6667 

1.8983 

2.2942 

t ‘ = 0.05 

t' = 0.10 

V B 0.15 

t' = 0.20 

1 1 A AP* 

■0 ■ B v.u;? 

.1 . f A jH 

t; • = u*jLU 

f /% n ^ 

u ■ = u .x;? 

j_ r ^ 

T> ■ s 

1.0053 

1.0217 

1.0510 

1.0969 

1.1280 

1.1664 

•1 1 .^ 

JU.itJUHV 

1.2739 

1.3510 

1.4534 

1-5955 

1 . 8 o <59 

£•1732 

2.9548 

1.0055 

1.0226 

1.0534 

1.1021 

1.1355 

1.1770 

•1 

1 . 295 T 

1.3830 

1.5016 

1.6715. 

1.9401 

2.1)231 

3*5554 

1.0057 

1.0235 

1.0557 

1.1069 

1.1423 

1.1867 

1.3155 

1.4121 

1.5456 

1.7409 

2.0589 

2.6513 

4 . 10 ll 2 

1.0059 

1.0243 

1-0577 

1.U13 

l,l 4 ^ 

1.1956 

n rtEC/T 
J-.Cr-JJfJ 

1.3337 

1.4388 

1.5860 

1 . 8 o 46 

2.1679 

2.8605 

4.60^ 

1.0053 

1.0217 

1.0511 

1.0972 

1.12^ 

1.1672 

1.2760 

1.3546 

1.4600 

1.6094 

l.Bim 

2.2614 

3.3327 

1.0056 

1.0228 

1.0539 

1.1032 

1.1373 

1.1799 

•1 

1.3033 

1.3961 

1.5259 

1.7212 

2.0524 

2.7440 

4.9304 

1.0058 

1.0238 

1.0566 

1.1093 

1.14^ 

1.1926 

1 

1.3312 
1.4392 
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1.8427 

2.2901 

3-3121 

6.9353 
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1.0248 

1.0592 

1.1152 

1.1547 

1.2052 

n r>*n *1 

J JLOr 

1.3594 

1.4831 

1.6677 

1.9704 

2.5455 

3.9357 

9.2308 - 
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TABLE II 

RATIO cm MOMESnS PCK COMERISSIBEE AND INCOMERESSIBLE ELO^S 


“1 



«o/«I 

■1 


t’ = 0.05 

BBa 

m/sm 

1^ 


1.0050 

1.0051 

1.0051 

1.0050 

1.0049 


1.0206 

1,0207 

1.0208 

1.0206 

1.0202 


1.0483 

1.0486 

1.0488 

1.0478 

1.0467 

.ItO 

1.0911 

1.0920 

1.0925 

1.0914 

1.0880 

.45 

1.1198 

1.1212 

1.1222 

1.1208 

1.1158 

•50 

1.1547 

1.1570 

1,1587 

1.1570 

1.1499 

•55 

1.1974 

1.2012 

1.2044 

1.2023 

1.1923 

.60 

1.2500 

1.2561(- 

1.2625 

1.2605 

1.2461 

.65 

1-3159 

1.3274 

1.3392 

1.3385 

1.3177 

.70 

1.4003 

1.4222 

1.4469 

1.4511 

1.4211 

•75 

1.5119 

1.5580 

1.6147 

1.6358 

1.5946 

.80 

1.6667 

1.7789 

1.9294 

2.0135 

1.97C7 

• 85 

1.8983 

2.2438 

2.7499 

3.1211 

3.1895 

• 90 

2.2942 

3.9315 

6.5741 

8.9709 

10.4271 


TABLE HI 

M07EMENT OB' CENTER OE ERESSURE AS FUNCTION OF STREAM 


MACH KUMEER AND THICKNESS COEHFICISIW 


“1 



c - 

c 

2a 






t’ = 0.10 

r — *"* — *' " 

t' = 0.15 

t' = 0.20 

|0.10 

1.0050 

-0.0001 

-0.0001 

-0.0002 

-0.0002 

.20 

1.0206 

-.0002 

-.0004 

-.0007 

-.0009 

.30 

! 1.0483 

-.0006 

-.0011 

- .0018 

-.0024 

.40 

1.0911 

-.0011 

-.0022 

-.0034 

-.0049 

.45 

1.1198 

- .0015 

-.0030 

- .0047 

-.0067 

.50 

1.1547 

- .0021 

-.0040 

-.0063 

-.0092 

.55 

1.1974 

-.0028 

-.0054 

-.0085 

-.0124 

.60 

1.2500 

-.0039 s 

-.0071 

-.0113 

-.0167 

.65 

1.3159 

- .0048 

-.0092 

- .0149 

-.0223 

.70 

1.4003 

-.0062 , 1 

-.0117 

-.0193 

-.0296 

.75 

1.5119 1 

- .0076 

-.0139 

-.0239 

-.0382 

.80 

1.6667 i 

-.0080 

-.0135 

-.0257 

-.0452 

•85 

1.8913 1 

-.0019 

.0005 

- .0123 

-.0381 

.90 

2.2942 ! 

.0427 

.0750 

.0624 1 

.0259 
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(a) Z-planei 
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(o) 


Pigur* 1.- Flow plane Z with directions of tangwt aM noirol on 

contour; plane * of affinely distorted profile; plane % of circle 
oonfonoalXy related to profile In e*-plano# 
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(a) t' = 0,05. COMMITTEE FOR AERONAUTICS 

Figure 2,- Ratio of lifts and ratio of moments in compressible 
and incompressible flows and movement of center of pressure 
in percent chord as functions of stream Mach number. Center 
of pressure movement rearweird with increasing stream Mach 
number. 
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Figure 2,- Continued, 
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Figure 2,- Continued, 
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Figure 2«- Concluded. 
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